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ABSTRACT 

We have previously reported on chemical abundance trends with evolutionary state in the 
globular cluster NGC 6397 discovered in analyses of spectra taken with FLAMES at the VLT. 
Here, we reinvestigate the FLAMES-UVES sample of 18 stars, ranging from just above the 
turnoff point (TOP) to the red giant branch below the bump. Inspired by new calibrations of 
the infrared flux method, we adopt a set of hotter temperature scales. Chemical abundances are 
determined for six elements (Li, Mg, Ca, Ti, Cr, and Fe). Signatures of cluster-internal pollution 
are identified and corrected for in the analysis of Mg. 

On the modified temperature scales, evolutionary trends in the abundances of Mg and Fe are 
found to be significant at the 2a and 3a levels, respectively. The detailed evolution of abundances 
for all six elements agrees with theoretical isochrones, calculated with effects of atomic diffusion 
and a weak to moderately strong efficiency of turbulent mixing. The age of these models is 
compatible with the external determination from the white dwarf cooling sequence. We find that 
the abundance analysis cannot be reconciled with the strong turbulent-mixing efficiency inferred 
elsewhere for halo field stars. A weak mixing efficiency reproduces observations best, indicating 
a diffusion-corrected primordial lithium abundance of loge(Li) = 2.57 ±0.10. At 1.2a, this value 
agrees well with WMAP-calibrated Big-Bang nucleosynthesis predictions. 

Subject headings: diffusion — globular clusters: individual (NGC 6397) — stars: abundances — stars: 
atmospheres — stars: evolution — stars: fundamental parameters — stars: Population II 

1. Introduction (2007, hereafter Paper I) and Lind et al. (2008, 

hereafter Paper II) - homogeneous, internally 

NGC 6397 is one of the best-studied metal-poor consistent analyses of stars at various evolution- 

([Fe/H] » -2.1) globular clusters of the Milky ary stages were p erforrned on spe ctra taken with 

Way. Prompted by recently determined external the mu iti_ bject spectrograph FLAMES at the 

constraints detailed below, we continue and up- VLT. In Paper I, the evolution of surface elemen- 

date our investigation of stars in this cluster. In a taJ abundances was traced from high-resolution 
series of articles - Korn et al. (2006), Korn et al. 
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FLAMES-UVES observations of 18 stars just 
past the turnoff-point (TOP), on the subgiant 
branch (SGB), at the base of the red giant branch 
(bRGB), and on the red giant branch below the 
bump (RGB). The elemental abundances of Li, 
Mg, Ca, Ti and Fe were found to display system- 
atic evolutionary trends incompatible with stan- 
dard stellar-evolution models, significant at the 
3a level in the cases of Mg and Fe. Comparing 
with stellar-evolution models including the effects 
of atomic diffusion as calculated by Richard et al. 
(2002, 2005), the trends were found compatible 
with those of theoretical isochrones for an age 
of <13.5Gyr, with turbulent mixing of moder- 
ate efficiency. These results were supported by 
the independent analysis in Paper II, where abun- 
dances were derived for Fe, Ti, Ca and Mg from 
medium-high resolution GIRAFFE observations 
of 116 stars along the subgiant branch in the same 
cluster. 

A number of recent studies have set new 
constraints on the temperature scale of metal- 
poor stars, cluster-internal pollution, the age of 
NGC 6397 and the efficiency of the turbulent mix- 
ing which moderates the effects of atomic diffu- 
sion. It is the aim of this article to reevaluate the 
abundance trends and their comparison to model 
predictions under these new constraints. 

On the observational side, the improved IRFM 
calibration of Casagrande et al. (2010) for dwarfs 
and subgiants produces ~ 100 K higher effective 
temperatures than the Ha analyses of Paper I. 
Similar alterations to the temperature scale were 
proposed by Bonifacio et al. (2007), as they would 
significantly reduce the strength of observed 
trends thus making the results of Paper I agree 
better with the analyses of Gratton et al. (2001) 
who found no abundance differences between TOP 
and bRGB stars. 

Na-0 anticorrelations as inferred by Carretta et al. 
(2009) from a very large set of red giants in 15 
globular clusters, indicate a largely bimodal dis- 
tribution of primordial (first-generation) and pol- 
luted (second-generation) stars. In NGC 6397, 
the primordial component in a sample of 16 RGBs 
with reliable oxygen abundances is found to en- 
compass roughly 25 % of all stars, similar to other 
globular clusters of the same study. For the same 
cluster, Lind et al. (2009b, 2011b) find a similar 
distribution of sodium abundances, tied to signif- 



icant anticorrelations with e.g. magnesium. 

The age of NGC 6397 has been determined 
in analyses of the white dwarf cooling sequence 
(hereafter WDCS) by Hansen et al. (2007) and 
Kowalski (2007). We adopt the latter's determina- 
tion as our WDCS age prior, 12.0 ± 0.5 Gyr (2a), 
to check whether the results of Paper I are modi- 
fied by this boundary condition on the theoretical 
modeling. 

The study of Gonzalez Hernandez et al. (2009) 
compares 3D-NLTE lithium abundances of main- 
sequence and subgiant stars at similar temper- 
ature in NGC 6397. They find a behavior on 
the upper subgiant branch in rough qualitative 
agreement with atomic-diffusion model predic- 
tions employing a high turbulent-mixing efficiency. 
Mclendcz et al. (2010) find good correlation with 
the mass-dependence of lithium depletion in field 
halo stars for diffusion models using the same 
(high) turbulent-mixing efficiency. Their deduced 
primordial lithium abundance, loge(Li) = 2.64, 
is in excellent agreement with predictions of the 
Big-Bang nucleosynthesis (BBN) lithium abun- 
dance using current WMAP-7 and nuclear reac- 
tion data and uncertainties, loge(Li) = 2.71 ±0.06 
(Cyburt et al. 2010). We evaluate whether mod- 
els adopting such high turbulent-mixing efficiency 
predict significant abundance variations in ele- 
ments other than lithium, compatible with obser- 
vations. 

In the present work, we revisit the stellar sam- 
ple and spectra of Paper I. We introduce three new 
temperature scales in § 2.2. We present our rede- 
termined elemental abundances on the new tem- 
perature scales in § 3, and our efforts to correct 
for effects of pollution in § 3.1. We compare these 
results to a grid of theoretical isochrones with ef- 
fects of atomic diffusion at various ages in § 3.3. 
We summarize and consider the implications for 
the cosmological lithium problem in § 4. 

2. Methods 

Photometric and spectroscopic analyses are 
performed as detailed in Paper I. We start from 
the spectroscopically determined, internally con- 
sistent stellar parameters determined there: MAFAGS- 
ODF ID LTE atmospheres (Grupp 2004) were 
used to iteratively derive T c g, log g, £, and [Fe/H] 
for each star. T e g was determined using Ha wing 
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Table 1 

Line list used in updated abundance determinations 



Species 


Line(s) [A] 


Treatment 


Li I 


A6707 


NLTE 


Nal 


AA5688, 8183 a , 8194 a 


NLTE 


Mgl 


A5528 


NLTE 


Cal 


AA6122, 6162, 6439 


NLTE 


Ti II 


AA5188, 5226 


LTE 


CrI 


A5206 


NLTE 


Fe II 


AA4923.9, 5018.4, 5169.0, 5316.6, 5234.6 


LTE 



a IR lines from supplemental spectra (Gratton et al. 2001; 
Lind ct al. 2009b). 



fitting employing Ali & Griem (1966) broadening 
theory, logg from the iron ionization equilibrium 
in NLTE (Korn ct al. 2003), and £ from the Fe 
line-strength trend. 

We deviate from this internally consistent ana- 
lysis by instead adopting absolute logg values 
derived from photometric V magnitudes, using 
the same parameters for reddening and distance 
modulus as in Paper I. We employ tempera- 
tures from the new temperature scales defined 
in § 2.2. Bolometric corrections are adopted from 
Worthey & Lee (2011). Stellar masses are taken 
from a 12.5 Gyr isochrone with atomic diffusion. 
Models without atomic diffusion would differ only 
by an absolute offset in logg. On the TOP and 
SGB, the expected ~ 0.25 dex underabundance in 
helium due to diffusion alters the mean molecular 
weight, corresponding to a decrease in spectro- 
scopic logg values (see Paper I). We account for 
this by decreasing logg values by 0.05 dex in our 
spectroscopic analysis of TOP and SGB stars. 
This assumption affects only the abundances of 
Ti and Fe (derived from ionic lines), decreasing 
both by a mere 0.01 dex. 

Chemical abundances are deduced from line-by- 
line solar-differential profile fits, thus cancelling 
the dependence on adopted logg/ values. Group 
average chemical abundances are not averages of 
individual stellar abundances - rather, they are 
derived from mean stellar parameters applied to 
group averaged spectra. This approach allows us 
to work on spectral lines as weak as 20 mA, as 



is the case for the ionic lines of titanium in the 
TOP stars. The line list employed in this updated 
analysis is given in Table 1. Abundance uncer- 
tainties on the differential scale are assumed to be 
0.05 dex, save for the cases of Mg (see § 3.1) and 
Fe (see Paper I, Table 3). 

2.1. NLTE and 3D modeling 

We make an effort to employ the most re- 
alistic line-formation theory available. NLTE 
corrections are applied for Mg I (Gchrcn et al. 
2004), Ca I (Mashonkina ct al. 2007) and Cr I 
(Bcrgemann & Cescutti 2010). Additionally, we 
apply NLTE corrections to LTE-dcrivcd abun- 
dances for Na I (Lind et al. 2011a) and Li I 
(Lind et al. 2009a). The majority species Fe II 
and Ti II are modeled in LTE, as they have 
been shown to exhibit only minor NLTE effects 
(Mashonkina et al. 2011; Bcrgemann 2011). 

We derive iron abundances from Fe II lines 
only, as the modeled strengths of Fe I lines de- 
pend critically on NLTE and its modeling assump- 
tions (specifically the inelastic H I collision factor 
Sh). As a diagnostic of this effect, we have also 
derived log g values from the iron ionization equi- 
librium on our new temperature scales. For this 
test, we employ NLTE corrections in the line for- 
mation producedure following Korn et al. (2003), 
with rather efficient hydrogen collisions (Sh = 3). 
These spectroscopically derived logg values over- 
all agree reasonably with those derived directly 
from photometry and distance. A significant dif- 
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ference is found on our hottest temperature scale, 
'200c', where log g values for the TOP stars are 
higher by 0.2 dex, placing them at logg = 4.2, 
which is just below the turnoff point in our theo- 
retical CMD (compare to Figure 1). 

3D effects on abundance trends in LTE were 
investigated in Paper I. Qualitative results implied 
stronger trends for Ti II and Fe II, but a weaker 
trend for Mg I. 

For chromium, we find similar NLTE correc- 
tions for the high-excitation triplet line Cr I A5206 
in all groups, ranging from 0.31 dex at the RGB 
to 0.37 dex at the SGB. Bonifacio et al. (2009) 
predict very large differential 3D effects at low 
metallicities. At [Fe/H] = —3.0, these correc- 
tions comparing giants to dwarfs are on the or- 
der of (5Aloge(Cr) = 0.4 — 0.5 dex for high exci- 
tation lines. Preliminary (3D) NLTE modeling 
implies individual corrections of similar magni- 
tude (M. Bcrgcmann, private comm.), rendering 
comparisons between dwarfs and giants uncertain. 
Thus, we expect that our results may change sig- 
nificantly under more sophisticated modeling. 

For lithium, we need not only consider differen- 
tial modeling effects, but must also determine ab- 
solute abundances to compare with the predicted 
primordial abundance. The 3D NLTE treatment 
by Sbordone et al. (2010), using an 8-level model 
atom, results in corrections A(3D — ID) = 0.03 dex 
for our TOP and SGB stars. In ID, these NLTE 
corrections (ANLTE) are very near zero. 

We have performed independent analyses fol- 
lowing Lind et al. (2009a) on metal-poor stars 
similar to our sample: for the metal-poor dwarfs 
HD 84937 and G64-12, and the subdwarf HD 
140283 representing well studied Spite plateau 
stars similar to our TOP and SGB stars, as well 
as an additional model representing a star evolved 
somewhat beyond our sample (at T e g = 5000 K, 
logg = 2.0, [Fe/H] = -2.0). In this analysis, 
ID models were constructed from MARCS at- 
mospheres (Gustafsson et al. 2008). Average 3D 
((3D)) models were constructed from sequences 
of snapshots from radiation-hydrodynamic con- 
vection simulations using the Stagger Code 
(Stein & Nordlund 1998), with averaging over 
the 3D cubes on surfaces of equal T500 (R- Col- 
let, private comm.). Results imply (3D) NLTE 
lithium abundances very similar to the corre- 
sponding ID NLTE case. On the Spite plateau, 



we find ANLTE = -0.05 dex in ID, with iden- 
tical or slightly higher abundances in (3D). For 
evolved stars, the correction instead becomes posi- 
tive. For the evolved model star investigated here, 
the NLTE abundance is lower in (3D) than in ID. 
The NLTE effects in this star should be greater 
than in our bRGB and RGB stars, as it has signif- 
icantly lower surface gravity. These results imply 
that hydrodynamical effects will not lead to sig- 
nificantly different results for our evolved stars, 
and that ID LTE analyses indeed underestimate 
their lithium abundances. 

Similarly small A(3D — ID) corrections result 
from these preliminary 3D and (3D) NLTE lithium 
abundance analyses. For this reason, we feel 
safe in applying ID NLTE corrections, but keep 
in mind the larger systematic uncertainty in the 
abundances of evolved stars. 

2.2. Scenarios 

We compare the effective temperatures de- 
duced from spectroscopic Ha analyses to those 
from three different photometric calibrations in 
Table 2. All photometric calibrations indicate a 
range AT c g(TOP — RGB) similar to or some- 
what more narrow than the spectroscopic analysis. 
Working on these temperature scales would gen- 
erally strengthen the previously identified trends. 

The Casagrande et al. (2010) calibrations for 
dwarfs and subgiants indicate higher effective tem- 
peratures for the TOP by 110 ± 30 K, but are not 
applicable to giants. Assuming the empirical cali- 
brations of Alonso et al. (1996) hold on the RGB, 
this would allow an increase in AT c ff (TOP — RGB) 
by 100 K. Note however that this is not the most 
likely interpretation of the hotter temperatures 
stemming from the Casagrande et al. (2010) cal- 
ibration. Rather, the difference is largely due to 
the calibration of the infrared bands, resulting in 
an absolute offset, rather than the inclusion of 
interfcrometric radii for dwarfs, which would not 
affect the giants and thus give a differential effect 
(L. Casagrande, private comm.). 

On the SGB, the effective temperatures are sim- 
ilarly higher by 115 ± 100 K, indicating no signifi- 
cant change in AT off (TOP -SGB). The errors are, 
however, fairly large, as the calibration sample is 
rather thin at this evolutionary stage (see Fig. 13 
of Casagrande et al. 2010). 
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Table 2 

Photometric effective temperatures 



Group 


Spectroscopic 




IRFM 1996/1999 a 






IRFM 2010 


b 


MARCS C 




Ha 


b — y 


v-y 


B-V 


V -I 


b — y 


B-V 


V -I 


b-y 


v-y 


TOP 


6254 


6229 


6213 


6160 


6133 


6374 


6373 


6327 


6288 


6281 


SGB 


5805 


5797 


5824 


5827 


5688 


5951 


6004 


5804 


5839 


5853 


bRGB 


5456 


5396 


5408 


5545 


5317 








5541 


5535 


RGB 


5130 


5121 


5132 


5254 


5063 








5225 


5226 


A(TOP - RGB) 


1124 


1108 


1082 


906 


1070 








1063 


1055 


A(TOP - SGB) 


449 


433 


390 


333 


446 


423 


369 


523 


449 


428 



Note. — Results have been averaged for each group. Photometric calibrations are evaluated for [Fe/H] = —2.0, 
on measurements listed in Paper I, Table 9. 

a Empirical calibrations for dwarfs (Alonso et al. 1999) and giants (Alonso ct al. 1996). We consider the SGB 
stars to belong to the dwarf calibration, and the bRGB stars to the giant one. 

b Empirical calibrations valid for dwarfs and subgiants (Casagrande ct al. 2010). 

c Synthetic photometry valid for dwarfs and giants (Onehag et al. 2009). 



The Ha analysis employs the broadening the- 
ory of Ali & Griem (1966). Direct comparison to 
that of Barklem et al. (2000) indicates only minor 
deviations, affecting AT cff (TOP - RGB) by 20 K, 
and increasing AT eff (TOP - SGB) by some 40 K. 

Starting from the spectroscopic temperature 
scale, we introduce three hypothetical tempera- 
ture scales. For each, photometric \ogg values are 
determined as described in § 2. In the most proba- 
ble scenario, we shift all effective temperatures by 
100 K, to match the expected absolute offset intro- 
duced in the Casagrande ct al. (2010) calibration. 
We refer to this temperature scale as '100s'. 

Additionally, we create two expanded tempera- 
ture scales. The first retains AT cff (TOP-SGB) in 
reasonable agreement with both photometry and 
spectroscopy, expanding the temperature scale lin- 
early, by a fixed amount for each stellar group as 
determined by its mean temperature. We expand 
AT eff (TOP - RGB) by 100 K, keeping tempera- 
tures fixed at the RGB, and refer to this temper- 
ature scale as T00e'. The increase in parameter 
AT eff (TOP - SGB), 40 K, is in agreement with 
the improved Ha broadening, and does not violate 
the photometric indicators. This scenario repre- 
sents a realistic but conservative case, relative to 



T00s'. The temperature scale '200e' is created by 
the same method, with AT cff (TOP - RGB) ex- 
panded by 200 K. This temperature scale, while 
poorly supported by our photometric and spec- 
troscopic analyses, is meant to represent the case 
where atmospheric parameters are tuned to flatten 
the abundance trends. 

The stellar parameters on each temperature 
scale are shown together with theoretical isochrones 
at three different assumed ages in Figure 1. 
The WDCS prefers an age of 12.0 ± 0.5 Gyr 
(Hansen ct al. 2007; Kowalski 2007), which the 
models effectively bracket. This age constraint is 
well fulfilled by all three temperature scales. Ad- 
ditionally, the spectroscopic values are compatible 
within their error bars. 

An additional constraint comes from the ob- 
served cluster CMD (see Paper I, Fig. 2). We 
expect our TOP stars to lie some 100 — 200 K 
cooler than the turnoff point (F. Grundahl, private 
comm.). This criterion disfavors the hot '200c' 
temperature scale as well as the spectroscopic de- 
termination - indicating a cluster too young or 
old, respectively - while T00s' and T00e' are both 
fully compatible. 

For traditional stellar models computed with- 
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6400 6000 5600 5200 

Effective temperature T cS [K] 

Fig. 1. — Comparison of predicted stellar parame- 
ters for theoretical T6.09 isochrones and the stellar 
samples. Symbols represent the various theoreti- 
cal temperature scales (see text). Shown as crosses 
are results with error bars from the spectroscopic 
analysis of Paper I. Note that the 13.5 Gyr model 
was computed at 0.1 dex higher metallicity, which 
produces a slight deviation on the RGB. The dif- 
ference is negligible near the TOP. 

out atomic diffusion, the lack of gravitational set- 
tling of helium results in a significantly older clus- 
ter. Then, only '200e' seems compatible with 
the WDCS age constraint, as well as the cluster 
CMD. From the cluster CMD constraint, '100s' 
and T00e' would both prefer ages greater than 
13.5 Gyr. 

2.3. Stellar models with diffusion effects 

Following Richard ct al. (2002, 2005), a sparse 
grid of isochrones has been computed with effects 
of atomic diffusion. Three different ages and six 
different efficiencies of turbulent mixing are con- 
sidered. At 11.5 and 12.5 Gyr, a metallicity of 
[Fc/H] = —2.1 was used. The 13.5 Gyr models 
were reused from Paper I, with [Fc/H] = —2.0. 
This has minor effects on e.g. the stellar mass 
at a given evolutionary stage, but is negligible 
for the evolution of surface abundances. Turbu- 
lent mixing is prescribed ad hoc with an efficiency 
parametrized by a reference temperature To deter- 
mining the behavior in both strength and depth 



(see section 3.2 of Richard et al. 2005, for more 
details). We use the shorthand Tx for models 
using log To = x, and refer to the magnitude of 
this parameter as the turbulent-mixing efficiency. 
The T5.95 and T6.25 models respectively repre- 
sent the lowest and highest turbulent-mixing effi- 
ciencies that seem compatible with the flatness of 
the Spite plateau of lithium. 

3. Results 

Mean stellar parameters and group-average 
abundances are presented in Table 3. 

3.1. Anticorrelations 

As noted in Paper I, the bRGB and RGB 
groups show some scatter in the magnesium abun- 
dance. Magnesium is known to anticorrelate with 
aluminum in globular clusters (see Gratton et al. 
2004), but no lines of the latter are detected in 
our spectra. We resort to using sodium as a 
proxy via its well-studied O-Na anticorrclation 
(e.g. Carretta ct al. 2009; Gratton ct al. 2004), 
and the related Li-Na anticorrelation which is 
known to be significant only for the most strongly 
polluted stars (Lind et al. 2009b). 

From a large sample of RGB stars, Carretta et al. 
(2009) have determined that some 25 ± 13 % 
of stars in NGC 6397 belong to the first gen- 
eration. The independent determination by 
Lind et al. (2011b) confirms this result. They 
identify for their sample of turnoff stars mean 
abundances for the first and second generations 
as loge(Na) « 3.65 and 4.05, respectively. For 
their combined sample of giants and diffusion- 
corrected dwarfs (+0.20 dex), the levels increase 
to loge(Na) = 3.84 ± 0.10 and 4.23 ± 0.13, re- 
spectively. As our sample size is insufficient for 
statistical analyses of the abundance trends, we 
rely on comparisons with the data and results 
of Carretta et al. (2009) and Lind et al. (2008, 
2009b, 2011b) for our population determinations. 
We perform the sodium abundance analysis on 
our original spectroscopic Ha temperature scale, 
as its good agreement with previous photometry 
makes it directly comparable with the cited stud- 
ies, giving us comparable abundance scales. This 
assumption is further supported by the agreement 
of stellar parameters for stars in common with the 
other studies. 
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Table 3 

Mean stellar parameters and abundances at the different temperature scales 



Group 


T e fi 


log g 


* [Fe/H] 


log e (Li) 


loge(Mg) b 


loge(Ca) 


loge(Ti) 


loge(Cr) 




[K] 


[cgs] 


NLTE 


NLTE 


NLTE 


NLTE 


LTE 


NLTE 




The 100 K shifted temperature scale TOOs' 


TOP 


6354 


4.01 


-2.23 


2.26 


5.80 


4.59 


3.05 


3.41 


SGB 


5905 


3.71 


-2.19 


2.39 


5.91 


4.58 


2.98 


3.46 


bRGB 


5556 


3.38 


-2.13 


1.46 


5.92 


4.63 


3.05 


3.49 


RGB 


5230 


2.58 


-2.07 


1.12 


5.98 


4.70 


3.10 


3.47 








The 100 K 


expanded temperature scale T00e' 






TOP 


6354 


4.01 


-2.23 


2.26 


5.80 


4.59 


3.05 


3.41 


SGB 


5865 


3.70 


-2.20 


2.35 


5.89 


4.55 


2.96 


3.42 


bRGB 


5486 


3.36 


-2.15 


1.41 


5.88 


4.58 


3.02 


3.41 


RGB 


5130 


2.54 


-2.10 


1.02 


5.96 


4.60 


3.05 


3.33 








The 200 K 


expanded temperature scale '200e' 






TOP 


6454 


4.03 


-2.21 


2.33 


5.84 


4.64 


3.09 


3.49 


SGB 


5925 


3.72 


-2.18 


2.41 


5.92 


4.59 


2.99 


3.48 


bRGB 


5526 


3.37 


-2.14 


1.43 


5.90 


4.61 


3.04 


3.46 


RGB 


5130 


2.54 


-2.10 


1.02 


5.96 


4.60 


3.05 


3.33 



Note. — Uncertainties in elemental abundances on the differential scale are assumed to 
be 0.05 dex. For iron, we use the statistical uncertainties of 0.04, 0.05, 0.04, and 0.03 dex for 
the TOP, SGB, bRGB, and RGB, respectively, as determined in Paper I. For magnesium, 
we estimate the uncertainty to be 0.065 dex for the TOP, SGB and bRGB groups, and 
0.05 dex for the RGB. 

a In the spectroscopic analysis, logg is corrected for the diffusion of helium by -0.05 dex 
for the TOP and SGB. 

b Magnesium has been corrected for pollution effects (see § 3.1). 
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For our sample of six RGB stars, we find a 
strong anticorrclation between magnesium and 
sodium, identifying two out of our six stars as 
belonging to the first generation in this cluster 
with loge(Na) w 3.9. Individual measurements of 
the very weak lithium line gives a relatively large 
scatter (<r w 0.10 dex), seemingly anti-correlated 
to sodium abundances at rather low significance. 

As indicated in Paper I, our sample of bRGB 
stars shows a bimodal distribution of magnesium 
abundances (A « 0.17 ±0.05). No anticorrelation 
is detected, as the- very weak sodium line seems 
to exhibit only mild scatter (er « 0.06 dex), with- 
out correlation to the former. Applying different 
effective temperatures from various photometric 
calibrations does not alter this result. Lithium 
shows significant scatter in line with expectations 
from an evolutionary scatter within the group (the 
cooler stars have more strongly diluted surface 
layers, with lower lithium abundances). Abso- 
lute sodium abundances (log£(Na) « 4.2) indicate 
that all five stars belong to the second generation. 

Within the SGB and TOP groups, no signifi- 
cant scatter is detected for magnesium. Sodium 
abundances identify both SGB stars as borderline 
second generation (loge(Na) » 4.0), while all five 
TOP stars clearly belong to the second generation 
(loge(Na) > 4.0). 

We account for pollution effects by attempt- 
ing to restore the magnesium abundances of the 
first-generation stars. We have directly identified 
first-generation stars in the RGB group, while for 
the other groups we must extrapolate from com- 
parisons with the aforementioned studies. We ap- 
ply upward corrections by 0.08, 0.07, 0.11, 0.11 dex 
to the group-averaged magnesium abundances for 
the RGB, bRGB, SGB, and TOP group, respec- 
tively. We propagate estimated uncertainties of 
the latter three corrections into those of the abun- 
dance determinations. 

If better abundance determinations for e.g. 
sodium or nitrogen were available, one could in 
principle reject polluted stars in the abundance 
analyses of sensitive elements. Unfortunately, it 
seems none of our TOP and SGB stars belong to 
the first generation, and thus may well represent 
a slightly decreased abundance in lithium. The 
magnitude of this effect may be gleaned from the 
observed difference amongst dwarfs in the large 
sample of Lind et al. (2009b), < 0.05 dex. Fu- 



ture studies could apply the identification method 
based on the nitrogen-sensitive photometric index 
c y when selecting their targets (see Lind et al. 
2011b, and references therein). 

3.2. Abundance trends 

On our preferred temperature scale, '100s', 
AT c ff is not altered, thus affecting abundance 
trends only by second-order effects. For mag- 
nesium, the corrections for pollution weaken the 
abundance trend to A log £ (RGB - TOP) = 
0.18 ± 0.08, which is significant on the 2a level. 
The trend in iron, meanwhile, is significant on the 
3ct level with Aloge(RGB — TOP) = 0.16 ±0.05. 
Calcium and titanium both show mildly signifi- 
cant trends, A logs (RGB - SGB) = 0.12 ± 0.07. 
Lithium displays a significant upturn at the SGB, 
Aloge(SGB-TOP) = 0.13±0.07, in line with the 
independent results of Lind et al. (2009b). The 
evolution of chromium is consistent with no vari- 
ation. 

Comparing to Paper I, magnesium abundances 
have been most significantly altered by correct- 
ing for cluster-internal pollution. This decreases 
A log e (RGB - TOP) by 0.03 dex, while increas- 
ing the uncertainty of this measure (from 0.07 to 
0.08 dex). The altered temperature scale has not 
affected abundance trends in neither magnesium 
nor lithium. Iron abundances are now based upon 
ionic lines, which are not temperature sensitive. 
The absolute abundances of iron and titanium are 
affected by the new, increased log<? values, with 
insignificant differential effects. The trend in cal- 
cium has strengthened by the combined differen- 
tial effects from T e g and logp by 0.04 dex. 

Expanding the temperature scale, i.e. in- 
creasing parameter AT c ff(TOP — RGB), implies 
higher abundances for the TOP stars. This leads 
to higher inferred abundances in temperature 
sensitive elements, which decreases parameter 
A log £ (RGB - TOP). This tends to flatten the 
observed abundance trends in Paper I, which all 
exhibited Alog£(RGB - TOP) > 0. 

On the expanded temperature scale T00e', 
trends in magnesium and iron are both significant 
on the 2(T level. Titanium and chromium exhibit 
similar but weakly significant trends of opposite 
sign. Calcium abundances are consistent with 
no variation. The upturn in lithium abundances 
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at the SGB remains, A log e (SGB - TOP) = 
0.10 ±0.07. 

On the hottest temperature scale, '200e', abun- 
dance trends by construction flatten significantly. 
The variation in iron however remains signifi- 
cant with A log e (RGB - TOP) = 0.11 ± 0.05. 
Chromium too deviates strongly, due singly to the 
behavior on the RGB. We note that while most 
signatures flatten, Aloge(TOP - SGB) of tita- 
nium actually increases somewhat on this temper- 
ature scale. This could be counteracted by arbi- 
trarily increasing logg values on the SGB, while 
inducing a complementary effect on iron. Flat- 
tening the iron trend further would similarly re- 
quire increasing log g values at the TOP, with the 
drawback of equally strengthening the trend of ti- 
tanium. Hence, this scenario results in nearly op- 
timally flattened abundance variations. 

3.3. Comparison with diffusion models 

In Figure 2, we compare observed abundance 
trends on the shifted temperature scale '100s' 
to stellar-structure predictions, as introduced in 
§ 2.3. Results on the expanded temperature scales 
are presented in Figure 3 and Figure 4. The 
dashed horizontal line represents the initial com- 
position of the models, adjusted to the observed 
abundance level. The precise choice of age, 11.5 
versus 12.5 Gyr, does not significantly influence 
the strength of the predicted evolutionary trends. 
Differences are larger when comparing to 13.5 Gyr 
models, in the sense that those predict stronger 
abundance trends for a given turbulent-mixing ef- 
ficiency. The influence of age is most significant 
on the cluster CMD - see § 2.2. From the dis- 
cussion there, we compare '200e' to the 11.5 Gyr 
isochroncs, and the other temperature scales to 
the 12.5 Gyr isochroncs, where more models are 
available. 

For the shifted temperature scale T00s', the 
models of least efficient turbulent mixing, T5.95- 
6.0, find excellent agreement with the observed 
trends in lithium, calcium, iron, chromium and 
titanium. For calcium, the predicted trend is 
somewhat stronger than observed, and for mag- 
nesium this is most apparent for the abundance 
at the SGB. Lithium, titanium and iron all ex- 
hibit trends favourable of weak efficiency. Taking 
into account the stronger upward adjustment re- 
quired to match the TOP and SGB abundances 



of lithium, the disagreement with the more un- 
certain RGB stars lessens when the T6.0 models 
are chosen. Hence, as in Paper I, this remains the 
optimal model, which gives the best overall agree- 
ment with observations. 

Models with moderately strong efficiency, 
T6.09 and T6.15, find mild disagreement with 
the predicted upturn in lithium abundances at 
the SGB, where instead a flat behavior is pre- 
dicted. Similarly, the evolutionary trend of iron 
becomes too weak. At less than 2a disagreement 
for all elements, these models can, however, not 
be excluded. 

Models with very strong turbulent-mixing effi- 
ciency, T6.25, disagree with the evolutionary trend 
of iron on the 2a level. Even worse, they dis- 
agree on the 3<7 level with the evolutionary trend 
of lithium, where a strong downturn is predicted 
at the SGB. This is because, with turbulent mix- 
ing this strong, gravitationally settled material is 
brought to sufficiently hot layers that lithium is 
destroyed. As the surface convection zone expands 
inward, it mixes with depleted material, thus di- 
luting the surface abundance further. Notably, 
this causes additional 0.15 dex of lithium destruc- 
tion, which naturally carries into the evolution 
along the RGB. 

Hence, models with turbulent-mixing efficien- 
cies in the range T5. 95-6. 15 agree well with obser- 
vations, with a preference for the weaker efficien- 
cies. A strong turbulent-mixing efficiency, T6.25, 
can be excluded for our sample of stars at this tem- 
perature scale from its predictions for both iron 
and lithium. 

On the expanded temperature scale T00e', 
the above results hold, with a possible prefer- 
ence for slightly stronger turbulent-mixing effi- 
ciency. The sole exception is chromium, where 
the TOP-bRGB range exhibits no variation, but 
Aloge(TOP - RGB) disagrees on the 2a level, 
with similar magnitude but opposite sign to pre- 
dictions. This result may however be significantly 
altered by (3D) NLTE modeling (see § 2.1). 

On the strongly expanded temperature scale 
'200c', trends by design flatten further, while not 
significantly altering the above conclusions. A 
weak efficiency of turbulent mixing is still com- 
patible with all observed trends, except for mag- 
nesium which clearly prefers somewhat higher ef- 
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Fig. 2. — Chemical abundances on the shifted temperature scale '100s', compared to predictions from 
12.5 Gyr isochrones. Dashed horizontal lines represent the initial composition. The stellar models include 
the effects of atomic diffusion with a free parameter for the turbulent-mixing efficiency, at five different 
values. We find that models in the range T5. 95-6. 15 all reproduce observations quite well. The optimal 
model is T6.0. 



ficiency, and chromium, discussed previously. 

4. Discussion 

We have updated the analysis of Paper I by 
deviating from its internally consistent spectro- 
scopic stellar parameters, and instead moving to- 
ward results inspired by the photometric calibra- 
tions by Casagrande et al. (2010). Further, we 
have corrected the abundances of magnesium on 
the grounds of two stellar generations in the glob- 
ular cluster, as identified by Carretta et al. (2009) 
and Lind et al. (2009b, 2011b). Only two of our 
stars may be easily identified as belonging to the 
first, pristine, generation, while the others seem to 
belong to a second generation with atmospheres 



to some degree polluted by more massive stars be- 
longing to the first. 

We have constructed three new temperature 
scales. Two of these find good agreement with 
our spectroscopic and photometric indicators: the 
shifted T00s' and the expanded T00c' tempera- 
ture scales. Both are compatible with the WDCS 
cluster age of 12.0 ± 0.5 Gyr (Hansen et al. 2007; 
Kowalski 2007). They exhibit evolutionary abun- 
dance trends, significant on the 2a level. On the 
shifted T00s' temperature scale, the abundance 
trend in iron is significant at 3(7. These results are 
compatible with predictions from stellar-structure 
models considering atomic diffusion moderated by 
turbulent mixing. Allowed models use turbulent 
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Fig. 3. — Chemical abundances on the expanded temperature scale '100e' (see Table 3), compared to pre- 
dictions from 12.5 Gyr isochroncs. Dashed horizontal lines represent the initial composition. 
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Fig. 4. — Chemical abundances on the hot expanded temperature scale '200e' (see Table 3), compared to 
predictions from 11.5 Gyr isochrones. Dashed horizontal lines represent the initial composition. 
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mixing of weak to moderate efficiency. The op- 
timal model is found to be the same as identi- 
fied in Paper I. The high-efficiency models pre- 
ferred elsewhere (Gonzalez Hernandez ct al. 2009; 
Melendez et al. 2010) produce additional destruc- 
tion of lithium, which leads to strong disagreement 
between our inferred abundances at the TOP and 
SGB. 

The third temperature scale, '200e', was con- 
structed under the criterion of optimally flatten- 
ing abundance trends. This is achieved, with the 
exceptions of iron and chromium at the 2a level 
(we dismiss the latter due to the large 3D cor- 
rections known to exist - see § 2.1). From the 
point of view of atomic diffusion however, we rec- 
ognize systematic albeit weak abundance trends 
in good qualitative accordance with predictions 
for all elements but chromium. Additionally, this 
temperature scale violates our spectroscopic effec- 
tive temperatures from Ha analyses, as well as 
all photometric calibrations considered. Finally, 
our magnesium abundances were corrected for 
pollution signatures deduced from sodium abun- 
dances. These were identified using the results 
of Lind et al. (2011b), where corrections for diffu- 
sion were applied. Neglecting this effect and in- 
stead calibrating on their large sample of dwarf- 
star sodium abundances results in a strengthened 
magnesium abundance trend significant at the 2a 
level. In this sense, it is not at all clear that there 
exists a temperature scale which does away with 
the need for atomic diffusion altogether. 

Finally, we acknowledge the differential ana- 
lysis of TOP and cool RGB stars in NGC 6397 
by Koch & McWilliam (2011). As our RGB stars 
have already undergone sufficient dilution of sur- 
face layers to restore all abundances but lithium 
to their original levels - signified by dashed lines 
in Figure 2 - one would not expect stronger effects 
to appear for their cooler stars. From their Fig. 4, 
with A = A(RGB - TOP), their absolute abun- 
dance variations in magnesium (A w 0.16 dex) and 
iron (A « 0.14 dex, based on Fe II) agree very well 
with ours, as well as model predictions. Their re- 
sults for titanium (A « 0.28 dex, based on Ti II) 
and calcium (A « —0.11 dex) deviate in alternate 
directions by 0.2 dex. 



TOP SGB bRGB RGB 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




6400 6000 5600 5200 
Effective temperature T eS [K] 

Fig. 5. — Observed abundance trends in lithium, 
shown for the shifted temperature scale T00s'. 
Lines represent the predicted surface abundance 
evolution at four different efficiencies of turbulent 
mixing. The horizontal dashed line represents the 
original lithium abundance of all four models, the 
dotted line the predicted primordial abundance. 
The shaded areas represent the respective uncer- 
tainties, and their overlap. 

4.1. The primordial lithium abundance 

In Figure 5, we show that the initial lithium 
abundance of the cluster as inferred from our opti- 
mal model (T6.0), loge(Li) init = 2.57±0.1, is com- 
patible with the predicted primordial abundance 
based on WMAP-7 data (Cyburt ct al. 2010, 
loge(Li)BBN = 2.71 ± 0.06). Adopting the mutual 
error bars, the values differ by 0.14±0.11 dex, or 
1.2a. If instead the T5.95 model were adopted, 
the initial abundance would further increase by 
0.04 dex to loge(Li) init = 2.61. Modeling in 
(3D) NLTE may further increase the abundance 
slightly, as would the assumption of some cluster- 
internal lithium destruction. Adopting a cluster 
age of 12.0 or 11.5 Gyr decreases the predicted sur- 
face depletion slightly, which reduces the inferred 
initial abundance by at most 0.01 dex. 

The model of high-efficiency model (T6.25) is 
not compatible with the lithium abundances on 
the TOP and SGB. The deep mixing causes ad- 
ditional destruction of lithium by 0.15 dex. When 
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calibrating solely on the RGB, this would result 
in the inferred initial abundance loge(Li) = 2.60. 
Note however that lithium abundances on the 
RGB are uncertain from the modeling perspec- 
tive as regards both line formation theory (see 
§ 2.1) and stellar evolution models (as indicated 
by a possible offset in T s or insufficient deple- 
tion of lithium along the subgiant branch, as seen 
in Figure 5, and likewise in Lind et al. (2009b)). 
These problems must be taken into account 
when considering the analysis of Mucciarelli et al. 
(2012). Indeed, these stars are less sensitive to 
the precise efficiency of atomic diffusion. But 
additional destruction of lithium occurs in models 
with strong turbulent mixing, an effect also affect- 
ing RGB stars, leading to underestimated initial 
abundances. We note in passing that their mean 
lithium abundance on the RGB for this cluster, 
loge(Li) = 1.09 ± 0.1, agrees with ours on the 
temperature scale T00s', loge(Li) = 1.12 ± 0.1. 

5. Conclusions 

We find agreement between the stellar parame- 
ters derived on two temperature scales compat- 
ible with the Casagrande et al. (2010) calibra- 
tion, and the stellar models with atomic diffu- 
sion of Richard et al. (2005) assuming an age for 
NGC 6397 of 12.0 ± 0.5 Gyr, as determined by 
Kowalski (2007). The spectroscopically deter- 
mined abundances of six chemical elements on 
these temperature scales exhibit significant varia- 
tions, in agreement with predictions from the stel- 
lar models if a weak efficiency of turbulent mix- 
ing is assumed. These results are rather robust as 
regards the choice of temperature scale. We can- 
not reconcile the detected abundance variations 
with predictions from models using the strong 
turbulent-mixing efficiency identified as optimal 
for Spite plateau stars in the field (Melendez et al. 
2010). 

On the preferred temperature scale, the op- 
timal model indicates an initial lithium abun- 
dance of loge(Li)i„it = 2.57 ± 0.1. This is com- 
patible with the predicted primordial abundance 
based on WMAP-7 data (Cyburt et al. 2010, 
loge(Li) B BN = 2.71 ±0.06). 

We note the WMAP-indcpcndcnt primordial 
lithium abundance of Stcigman (2010), loge(Li) = 
2.65 ± 0.06, in excellent agreement with our re- 



sults. This value is computed from the baryon 
density compatible with deuterium abundances in 
high-rcdshift QSO absorption line systems. Addi- 
tionally, one should expect some reduction from 
the primordial abundance due to mixing through 
population III stars (a mild version of the scenario 
presented by Piau et al. 2006). Finally, the first 
and second generation stars in this cluster may 
exhibit slightly different lithium abundances. Pre- 
vious studies (Lind et al. 2009b) identify lithium 
abundances higher in first than second generation 
dwarfs by < 0.05 dex. We note that a differential 
analysis comparing our TOP stars to the field star 
HD 84937 implies equal lithium abundances. Fu- 
ture studies could focus on analyses of first genera- 
tion stars, to remove this minor possible bias when 
constraining the true initial lithium abundances of 
cluster stars. 

In view of these uncertainties on either side, 
one may call into question the significance of the 
remaining cosmological lithium discrepancy. We 
can, however, be more certain about the reality of 
atomic diffusion in stars on the traditional Spite 
plateau, as shown by theoretical and observational 
investigations alike. What is still missing is an un- 
derstanding of the physical processes that give rise 
to the mixing required to moderate atomic diffu- 
sion. We call on stellar-structure experts to inten- 
sify investigations into possible (hydrodynamic) 
effects. 
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